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Abstract 

We transform an inverse scattering problem to be an interior transmission problem. We find an 
inverse uniqueness on the scatterer with a knowledge of a fixed interior transmission eigenvalue. By 
examining the solution in a series of spherical harmonics at far fields, we can decide the perturbation 
uniquely for the radially symmetric perturbations. 
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1 Introduction 

In this work we study the inverse acoustic scattering problem of recovering the index of refraction in an 
inhomogeneous domain. To determine the inhomogeneity, we send a wave field into the domain. The 
propagation of the detecting wave field will be perturbed when hinging on the inhomogeneity defined by 
the index of refraction that in turn produces a scattered wave field. The inverse problem is to determine 
the index of refraction by the measurement of the scattered wave field. The study of inverse scattering 
problem is the core in various disciplines of science and technology such as sonar and radar, geophysical 
sciences, medical imaging, remote sensing, and non-destructive testing in instrument manufacturing. 

In this paper, we take the incident wave field to be the time harmonic acoustic plane wave of the 
form 

u\x) := e ikx ' d , 

k £ R + , x £ R 3 , and d £ S 2 is the impinging direction. The inhomogeneity is defined by the index of 
refraction n £ C 2 (R 3 ); n{x) = n(|a;|) > 0 and n(|x|) — 1 y 0, for x £ ft, a ball of radius R in R 3 . The 
wave propagation is governed by the following equation. 

{ A u(x) + k 2 n(\x\)u(x) = 0, x £ R 3 ; 

u(x) = u l (x) + u s (x), x £ R 3 \ Q; (1.1) 

lim ixi-s-oo MlAiir “ ikuS ( x )I = °- 

Particularly, we have the following asymptotic expansion on the scattered wave field [HE]. 

£ ik\x\ ^ 

u s {x) = -j-j-UooOz; d, k ) + 0{— 3 -), |x| —> oo, (1.2) 

M |z|2 

which holds uniformly for all x := |^t, x £ R, and itoo(£;d, k) is known as the scattering amplitude in 
the literature pi nil ei Hill]- We state the following inverse uniqueness result in this paper. 

Theorem 1.1. Let u^^x; d, k), j = 1,2, be the scattering amplitude parametered by the non-constant 
index of refraction rP £ C 2 (R 3 ) in m- If u] X) {x\ d , k ) = w^ 0 {x\ d , k) for all x £ § 2 with a fixed d £ § 2 
and a fixed k > 1 , then n 1 = n 2 . 
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Let u SJ ( x ) be the scattered wave field parametered by n J . From Rellich’s lemma in scattering theory 

eh eh eh, the Theorem o assumption and HU imply 

u sl (x ) = u s2 (x ), i€l 3 \ 17. (1.3) 

Most importantly, (HU), (HU, and HU imply the following interior transmission problem. Let us set 

w(x) ■.= u sl (x), v(x) := u s2 {x ), 


and thus we have the following system of equations. 


Aw + k 2 n}w = 0, in 17; 
Av + k 2 n 2 v = 0, in 17; 
w = v, on 917; 


dw _ dv 

dv ~ dv' 


on 917, 


(1.4) 


where v is the unit outer normal. The equation HU is called the homogeneous interior transmission 
eigenvalue problem El IH IH EH EH EH EH EH EH EH EH EH EE] ■ The problem (|1.4I) occurs naturally 
when one considers the scattering of the plane waves by certain inhomogeneity inside the domain 17, 
defined by an index of refraction in many models. 

Now we expand the solution (w, v) of HU individually in two series of spherical harmonics by 
Rellich’s lemma [14] p. 32]: 


1 oo m=l 

w(x\ k) = - Y. Y. ai, rn ai(r\k)Y l m {x ); 

l —0 m=—l 
oo m=l 

v{x;k) = ~ Y Y 

l —0 m=—l 


(1.5) 

( 1 . 6 ) 


where r := |a;|; x = (0,^) € S 2 , k £ C. The summations converge uniformly and absolutely on suitable 
compact subsets in |r| > i?o, with some sufficiently large Rq > R. The expansion holds uniquely in the 
exterior domain of Helmholtz equation. The spherical harmonics 


{Yr(6,<p)h, m :=U 


2 4t 1 SI+|:!!l ^ l(cos ^ } -- m= ~ i > ■ ■ ■ ,i ' 1 = °’ ^ 2> • • ■ ’ 


(1.7) 


form a complete orthonormal system in L 2 (§ 2 ), in which 

Pn(t ) := (1 - *T /2 ^r^, m = 0,1,..., n, 

where the Legendre polynomials P n , n = 0,1,..., form a complete orthogonal system in L 2 [— 1,1]. 
By the analytic continuation of Helmholtz equation, the expansions HU and (11.61) converge up to the 
boundary dfl, that is, |a;| = R. 

The orthogonality of the spherical harmonics EH p. 227] implies that the functions in the form, 

( a itm (x-,k) := 

{ b liTn (x;k) := /3i, m bi( ]? fc) T; m (z), 

satisfy the first two equations in HU independently for each (l,m) for |x| > R. To fulfill the boundary 
condition of HU, we look for any k £ C such that 


( ai,mai(r\k)\ r = R = /3i, m bi(r; k)\ r=R ; 

l a l , m d r ^\ r=R = Pl, m d r b -±Y 1 \ r = R . 


(1.9) 
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By linear algebra, the existences of ap m and (3p m in (11.91) are equivalent to finding the zeros of the 
following functional determinant: 


Di(k\ R) := det 


ai{r;k)\ r=R bi(r\k)\ r=R 

( ai(r;k) / r bi(r;k) ~i / 1 

V y J h —H \ 7" s |r—^ 


that is, 


Di ( k . = ai(R; k)b[(R; k) - a'^R; k)bi(R; k) ^ 1Q ^ 

R 

Due to the radially symmetric assumption on m, the Fourier coefficients ai(r;k) and bi(r;k) solve 
the following system of ODE for all ( G No: 

( a"(r; k) + (k 2 n x (rx) - l< ' l ^ 1 ' > )ai (r; k) = 0, 0 < r < oo; 

\ 6"(r; k) + (k 2 n 2 (rx) — ^^-)bi(r; k) = 0, 0 < r < oo; (1-11) 

{ Di(k-,R)= 0. 


By the uniqueness of the Sommerfeld radiation condition of (HU), we deduce that 


&l,m — — 1 ■ 


( 1 . 12 ) 


We also set the initial conditions of aj(r; k) and 6;(r; k) at r = 0 to be the following conditions. 

- ji{kr)} = 0; (1.13) 

i—>-o r 

lim{ — ’ — ji(hr)} = 0. (1.14) 

i—>-o r 

The behavior of the Bessel function ji(kr) near r = 0 is found in [211 p.437]. We refer the initial 
condition (11.131) and (11.1411 to [231. Considering (11.131) and (11.141) , we deduce (11.111) to be the following 
ODE. 

a"(r;fc) + {k 2 n 1 (r) - i ^ 1) )a;(r;fc) = 0, 0 < r < oo; 

< b'/(r;k) + (fc 2 n 2 (r) - l -^ L )bi{v,k) = 0, 0 < r < oo; ^ 

Di(k; 0) = 0; 

. D l (k;R)=0. 

This is the reduced inverse problem ED in radially symmetric form. Because v)}i,m is complete 

in L 2 (§ 2 ), one element of the basis can be replaced by an element of fractional order of ( l,m ). Hence, 
the property m holds for l > 0, m < |/| without loss of generality. 


2 Asymptotic Expansions and Cartwright-Levinson Theory 

To estimate the asymptotic behaviors of the solution ai(r;k), bi(r;k), and then Di(k]R), we consider 
the following Liouville transformation: 

zli^lk) := [n 1 (r)]*ai(r;k), where ^ := [ [^(p^dp; (2.1) 

Jo 

z?(f-,k) := [n 2 (r)]*bi(r;k), where £ 2 := [ [n 2 (p)]^dp. (2.2) 

Jo 


Let us set r 

B J ■■= [ [n° (p)] 5 dp, j = 1,2. 

Jo 

If k is an interior transmission eigenvalue, then 

r $}" + [ k 2 - p'it' Yzj = 0, 0 < & < B\ j = 1, 2; 
1 A(fc;0) = 0; Di(k; R) = 0, 


(2.3) 
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in which 


P>(?) ■■= 


[n J ]"(r) 5 {{n 3 ]'{r)} 2 1(1 + 1) 


(2.4) 


(2.5) 


4 \ni(r)} 2 16 [n-l(r)] 3 r 2 r+(r) 

Here = B ,J if and only if r = R. Let 

= 5 {H'(r)} 2 Z(Z + 1) _ 1(1 + 1) 

J ' 4[nl(r)] 2 16 [n-?(r)] 3 r 2 nl(r) [£?] 2 ' 

Thus, (12.31) and (12.41) imply that 

f [+}" + [ k 2 _ gjtf) _ |+i)] z i = 0, 0 < < IB: 

\ A(fc;0) =0; Di(k; R) = 0. 

Let us drop the superscripts for simplicity. When l = 0, the estimates for the solution zo(i’) are 
classic and can be found in [251: 


sinfc£ cos kf... sin k£, ... 1 _,.exp[|Srfc|£]. 

2 °^; fc ) = nr - + 1^^) +P(°) - 2 Q (6] + ° ( fc4 )> 


where Q(£) := f 0 p(s)ds and the error term has an improvement [55J p. 17] by 


cos , 
8fc 4 


{V(£) -p'(0) - [p(£) +p(0)]Q(0 - J p 2 (s)ds + ^Q 3 (C)} + o( 


exp 

1*1 4 } ‘ 


Similarly, 


. . sinfc£... cosk£ r ... _,exp[|9ffe|£]. 

-o(£; fc ) = cosfc C + -^r-QCO + -^2 ^[p(C) -p(o) - (0] + °(-)» 


( 2 . 6 ) 


(2.7) 


( 2 . 8 ) 


in which the boundary behavior of p(£) plays a role in determining the inverse spectral uniqueness on 
the scatterer, and thus the C 2 - assumption on the index of refraction is necessary. 

For l > — i, we apply the much more generalized results from M: Let zj(£, k) be the solution of 


-*f (0 + 

lim^-o jrrr < oo, 


(2.9) 


in which the function g(£) is assumed to be real-valued and square-integrable. We note that the initial 
condition (11.131) and (11.141) imply the regularization condition at £ = 0 in (12.91) . The following asymptotics 
hold [7] Lemma 3p. 855]. 

. ... 7 . sin Jvlog(l + |fc|) 

I zi&k) - < - 6V 1 exp{|^fe|g}; 


W(C,k) - cosfcfl < Kl0g ^ + ^ exp{|£yfe|^}, where K = /v(||g||). 

1*1 

This explains the behaviors of solutions z;(£; k) for all l > — 

For the special case that q(£) = 0, we are actually considering the Bessel’s equation: 

< + [ k 2 - = 0 . ( 2 . 10 ) 

The solutions of (12.101) essentially are Bessel’s functions with a basis of two independent elements. 

The variation of parameters formula leads to the following pair of integral equations connecting 

zi(£,k), ui(£,k): 

zi(€,k) =ut(£,k)- f G(£, t, k)q(t)zi(t, k, q)dt, 

Jo 
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where 


G(£, t, k) = 

in which $(z,w) = <j>\{z)<t> 2 {w) — </>i(w)</> 2 ( 2 ) is defined as in [5] p. 6] that satisfies <J>(w,w) = 0 and 

7TT 1 

Mx) = (—)^Y l+h ( x y, 
fa (z) = (y)^ +i (A), 

where J v {z) is the Bessel function of the first kind and Y v (x) is the Bessel function of the second kind. 
Similarly, we have the following integral equation. 

z i(& k ) = u'li&k) ~ [ H(£,t,k)q{t)zi(t,k,q)dt , 

Jo 

H(t;, t, k ) = fc _1 \ i(kt, k£), 

in which T^w) = — <j>\(w)<l> 2 {z) with is defined as in [5] p. 6] which satisfies ^(uqw) = 1. A 

pair of solutions of (12.101) , are given by 



ui(£; k) = k l (j) i(fc£); 

(2.11) 


u 2 (^,k) = k~^MkO- 

(2.12) 

Moreover, we recall that 

jl ( z ) = 

(2.13) 

which we refer to [231 P- 

437], wherein we find that J u (z) and Y u (z) are holomorphic functions of z and 

entire functions of the order v when z is fixed. In this paper, the complex analysis 

V. 

For ^ > 0 and 9f?fc > 0, there is a constant C such that 

is focused at 


Mf; *) - Sinl jf + 7 1 1 < c|t|-< ,+1 > exp []*) KI ; 

(2.14) 


K(£;*)- cos{ lr' f) liC|* : r‘ exp ]^ } . 

(2.15) 

Moreover, 

I*i(6 k) - Ul (ik) 1 < C{—) l+1 exp{|3fc|£}£(£; fc); 

1 + \k£\ 

(2.16) 


wm k) - «£(€; k) 1 < exp{|S/c| £}E(b k), 

(2.17) 

where 

E({;t)=exp{^ 1+ y (| <B> 1. 

(2.18) 


We refer these estimates to [BJ Lemma 2.4, Lemma3.2]. These estimates give the fundamental asymp¬ 
totic behaviors of o;(r;fc), a[(r;k), bi(r\k ) and 6j(r;fc) and, ultimately, the behavior of Di(k',R). Most 
important of all, the estimates from (12.141) to (12.171) show that they are entire functions of order one 
and of type £. We refer details to [DEED! 02 EES Hang. Accordingly their zero sets are described by 
Cartwright’s theory ummmm- In particular, this leads to Weyl’s type of asymptotics for the zeros 
of Di(k) which we describe as follows. 

Definition 2.1. Let f(z) be an integral function of order p, and let N(f,a,/3,r) denote the number of 
the zeros of f(z) inside the angle [a,/3] and |^| < r. We define the density function of the zero set as 

N(f,a,p,r) 


A f(a, /3) := lim 


fP 


and 


A f{@) ■= A/(a 0 ,/3), 

with some fixed ao £ E such that E is at most a countable set 
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Theorem 2.2. The functional determinant Di(k; R) is of order one and of type B 1 + B 2 . In particular, 


ADi(k-,R)(—e>e) — 


B 1 + B 2 

IT 


Similarly, 


a / . B^ . . . B 1 

^ai(k-,R)\—C,£) — -! A a /( fc;i j)( —C, e ) ~ -i 

7 r 1 7r 

td2 td2 

e, e) = — ,A b 'i k . R ){-e,e) = —. 

7T 1 7T 


Proof. To sketch an idea of the proof, we note that the growth rate of Di(k; R) comes either from 
ai(R\ k)b[(R; k) or a[(R; k)bi(R; k) by considering (11.101) . Either of them have growth rate B 1 + B 2 if we 
examine the estimates from (12.141) to (12.171) . These two terms do not cancel each other if n 1 ^ n 2 . We 
skip the details and refer it to nnnmuMS]. 

□ 


In particular, the spectrum of dm is not empty and discrete in C. 


3 Proof of Theorem 11.11 


Let k be an eigenvalue of El- From the Rellich’s expansion E 3 and (II.CD . (11.91) and then (11.151) hold 
for all l > 0. To prove Theorem 11.11 we consider an inverse Sturm-Liouville problem to (11.151) with a 
common eigenvalue holding for all l > 0. 

We recall that R 7 and zj (£ J ; k ) are the quantities parametered by the index of refraction n - 7 , j = 1,2, 
as in (12.21) and (12.(11) . Therefore, the estimates (12.141) . (12.151) . (12.161) . (12.171) . and (12.181) hold for all Z > 0 
and for both parities of indices of refraction. We have 






(3.1) 


and 

no 

\zi(B^-k)- Ul (B 2 -k) I < C( 1 + ) l+1 eMmmE(B>-,k), (3.2) 

in which exp{|Q : fc|i? J } are bounded in a strip S containing the real axis, and E(B^,k), B J £ R + , is 
decreasing in k. Firstly we apply ED to obtain 


ui(B^-k) 


sin {kB> - l f} 
k l+1 


[l + 0( 


1 

\kBi\ 


)], l £ N 0 , 


(3.3) 


in which k is not in the zero set of sin{fcR- 7 — ^ } and the function is bounded near the real axis. Moreover, 
we deduce the following formula from (12.121) and (12.131) . 


7T Jc /~? ^ 1 

k l+1 ui(B j ; k) = (—— J l+ i(kBi) = kB’JuiklP). (3.4) 

2 2 

We note that a function of the form zji(z) is called a Riccati-Bessel function [24] . 

Furthermore, 

\k l+1 Ui(B l \ k) — k l+1 ui(B 2 -k)\ < k) — k)\ 

+\k l+ 1 zi(B 1 -k)-k l+ 1 zf(B 2 -k)\ 

+| k l+1 zf(B 2 ; k) - k l+1 ui{B 2 ; fc)|. (3.5) 

From (1321) . 

\k l+1 zi(Bi-,k)-k l+1 Ul {Bi-k)\ < ) l+1 e^ kBi 'E(Bi-,k) = o(l), as l -»■ oo, (3.6) 
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in which the estimate of E(B • J ; k) is given in (12.181) . 

Secondly, an eigenvalue k of d is surely an eigenvalue of (11.151) and (11.91) . By Liouville’s transfor¬ 
mation urn m with n 1 (i?) = n 2 {R), we deduce for this k that 


k) = zf(B 2 -, k), VI G R + . 
Assuming k > 1 by Theorem 11.11 assumption. (13.51) . and (13.71) imply 

|fc i+ 1 u;(i? 1 ; k) — k l+1 ui(B 2 -,k)\ —> 0, as l — > oo, 


that is, 


| (kB 1 )ji(kB 1 ) — (kB 1 )ji(kB 2 ) | — > 0, as l — » oo. 
Lemma 3.1. If H3.8\ ) holds, then B 1 = B 2 . 

Proof. From the Wronskian Identity [23J p. 439, (10.1.32)], we have 

jn+\{z)y n -\{z) - j n -i(z)y n+ i(z) = (2n + l)z” 3 


(3.7) 


(3.8) 


Because y n (z) = (—l) n+1 j- n -i(z), neZ, we obtain that 

(-!)" jn+l{z)j-n(z) ~~ {-l) n jn-l(z)j-n- 2 (z) = (2n+l)z~ 3 . 


(3.9) 


Moreover, 


Y Az) 


J v {z) COs{^7r} — J_v(z) 
sin{i/7r} 


Let us restrict v £ Z. Thus, J v (z ) = (—1 ) V J- V (z), and j l/ _i(z) = (— l) v j_ v _ i ( 2 ) accordingly by (12.131) . 
Hence, (13.81) is deduced to be 


|(fcB 1 ). 7 I (fe.B 1 ) - (kB 1 )ji(kB 2 )\ -»■ 0, as l 


±00 in — + Z. 


(3.10) 


In general, j v {z) is entire in ^ and | J v (f$tz)\ < 1 for v > 0 [Ml (9.1.60)]. There are convergent subsequences 
for {(kB 1 )ji(kB 1 )} l& i +z and {(fcH 2 )jj(A:i ? 2 )} 26 i +z with kB 1 and fcH 2 hxed. Hence, the convergence 
holds for l £ M + , then for l £ R. Then, applying (13.91) for both parities imply that 


(-1 ) n Jn+1 (kB l )j_ n (kB l ) - (-1 ) n j n _ 1 {kB i )j_ n _ 2 {kB i ) = (2n + l)(kB i )~ 3 , i = 1,2, 

that is, | 21 — 2 i fc T ) ? t 1 1 —> 0, as m —> 00 for some subsequence {7 m } G N. This is not possible unless 

B 1 = B 2 . 

□ 


Applying this lemma, now we want to show that z^B 1 ; k) = £q(.B, k) in C. In general, (13.21) implies 
that 

zi(B 1 -,k) = ui(B 1 -,k) + 0[( I -^ r ) l+1 eMmB 1 \}E(B 1 -,k)}, j = 1,2. (3.11) 

We apply (13.31) to (13.111) . and obtain 


zj(B 1 -k) = 


sin {kB 1 - l f} 
k l +i 


[l + 0(-)], i = 1,2. 


Thus, we deduce the following asymptotic behavior of the quotient. 


zl{B\k) 

zf(Bi-,k) 


i + Q(i) 

1 + o(i) 


— 1 H - 0( —), & G C, 
k 


(3.12) 


(3.13) 


sinj* _ Izl \ _ 

outside some neighborhoods of — fci+1 2 . Particularly, for l = 0, we apply the asymptotics (12.61) 

to (13.131) . Let { 7 ,;}^! be the simple zeros of sln H B 1 and {r j }^T 1 be some sequences of neighborhoods 
containing { 7 j}‘-fL 1 . We require that |Tj| —> 0, as j — » 00 . 
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For k {r j }°l 1 , we note from (13.1311 that 


,. zp( Bl ;k) 

k-toozjftB^k) 


For k G {Fj}^, we consider the following quotient. 


(3.14) 


lim lim 


(B x ;k) 


= lim 


sin{kB 1 y/k 


j^oo fc-> 7j z^B 1 ; k) j-koo, |r,|/o q 1T1 fc , /; 


kj sin{ kB 1 }/k 


From m an d (|2.7p , we have the following power series. 


zb( Bl > k ) = l _ CQt kB 1 q -j.) + [p{ b1 ) +P(0) - \Q 2 ( Bl )\ 


sin {kB 1 } 
k 


2 k 


4 k 2 




Moreover, by complex analysis, we have 

z bi Bl ;k) Res {zl(B l -k)l sm{k k B } ; 7 j} 


sin-ffcB 1 } 
k 


i k - 7 j) 


+ higher order terms , for k near 7 j. 


Now (13.161) implies that 


Res{^o(R 1 ;fc)/ Sm ^~ g ^ ; Ty} 

= Res{l - COt ^ Q\ Bl ) - gF[ q'(B 1 )} 2 + ... ;t,-} 


1 


x + 0 ( 4 )- 


Note that Res{ sin |^, B i| ; 7 j} = (—1 Y/B 1 . Hence, 

sialkB 1 } 


Re 8(2 J(B>; W ^i7^; 7j } = -%£l + o(±). 


k ’ IJJ 2S 1 7j 
Let us examine Q 1 (B 1 ). By integration by parts, we deduce that 

, -B 1 1 l" 

1 n 


Ti 


Q 1 (B 1 )= [ 

Jo 


4/ 


5 n 1 ' , 3 

*=16 


B l' 

n 


16 „i 3 


Q ds > 0 , 
13 — ’ 


(3.15) 


(3.16) 


(3.17) 


in which we emphasize that n 1 ^ 1 = B 1 ) = n 1 (r = R) = 0 by the fact that n 1 G C 2 (R 3 ) and is 
constant outside 0. Here, Q 1 (B 1 ) = 0 if and only if n 1 = 1. The same analysis holds for the index n 2 as 
well. Assuming n 1 and n 2 are non-constant as the assumption of Theorem 11.11 (13.151) and (13.171) , with 
non-vanishing second coefficient Q 1 {B 1 ) and Q 2 (B X ), imply that 


Hence, (13.141) and (13.181) imply that 


Bm Bm iflflifi = 1 . 

i->oo Zq (£>,£;) 


lim = 1. 


fc->oo z^B 1 ; k) 


(3.18) 


(3.19) 


This implies that -zq-frnTx has only finite number of irreducible zeros, denoted as {zl, z\, ■ • •, zlA, or 

Zq{Jd ]K) 


2 2 2 
Zi,Z 2 ,...,Z A 

zeros and poles are equal to some M G Nq. Let 


poles, {zf, z 2 ,..., zj^}, in 0 i + M, in which we deduce from (13.191) that the numbers of the irreducible 


F(k) := 


z o{B x 'i k)' 
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Therefore, 


(3.20) 


= (fc - z\)(k -zp---(k- zlf) 
(k~ zf)(k- z$)---(k- z 2 M )' 


We note that whenever k is an interior transmission eigenvalue of (11.151) . then it satisfies = 

4(B\k) by (HU), (irm (HU , and (22). That is 


F{k) = 1. 


Theorem 12.21 implies that 

B 1 + B 2 B 1 B 1 +B 2 B 2 

- > —, - > —, 

7r 7T 7T 7r 


and then there is a higher density of zeros of interior transmission eigenvalues than the Dirichlet eigen- 
values or Neumann eigenvalues of (12.91) . Hence, we deduce from the Fundamental Theorem of Algebra 
to (13.201) that 

F(Je) = 1. 


In particular, Zq (£; k) and Zg(£;fc) have the same Dirichlet eigenvalues. 

Similarly, we can prove that [zo]'(£; k) and [^o] , (^5 &) have the same Neumann eigenvalues by consid¬ 
ering (12.81) and 


l4]'^-k) 

k) 


— 1 + 0( —), k G C, 


outside some neighborhoods of the zeros of cosjfci? 1 }. If n 1 and n 2 have the same set of Dirichet and Neu¬ 
mann eigenvalues, then the inverse uniqueness result of the Bessel operator [BJ Theorem 1.2, Theorem 1.3] 
implies that n 1 = n 2 . This proves Theorem ll.il 


□ 


Remark 3.2. It is believed that a stepwise potential function in general can not be recovered by one 
discrete spectrum. See the Livshits’ example on potentials in this class m- In this paper, a C 2 -index of 
refraction defines a series of ODE from far-fields to the origin. 
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